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ABSTRACT 

We investigate the stellar populations of Lya emitters (LAEs) at z = 5.7 and 6.6 in a 0.65 deg^ sky of 
the Suharn/ XMM-Newton Deep Survey (SXDS) Field, using deep images taken with Subaru/ Suprime- 
Cam, UKIRT/WFCAM, and Spitzer/IRAC. We produce stacked multiband images at each redshift 
from 165 {z = 5.7) and 91 {z = 6.6) IRAC-undctectcd objects, to derive typical spectral energy 
distributions (SEDs) of 2; '--^ 6 - 7 LAEs for the first time. The stacked LAEs have as blue UV 
continua as the HST/WFC3 z-dropout galaxies of similar Muvj with a spectral slope /3 ^ —3, but 
at the same time they have red UV-to-optical colors with detection in the 3.6/xm band. Using SED 
fitting we find that the stacked LAEs have low stellar masses of ^ (3 — 10) x IO^Mq, very young ages 
of ~ 1 — 3 Myr, negligible dust extinction, and strong nebular emission from the ionized interstellar 
medium, although the z ~ 6.6 object is fitted similarly well with high-mass models without nebular 
emission; inclusion of nebular emission reproduces the red UV-to-optical colors while keeping the UV 
colors sufficiently blue. We infer that typical LAEs at 2: ~ 6 — 7 are building blocks of galaxies seen 
at lower redshifts. We find a tentative decrease in the Lya escape fraction from z = 5.7 to 6.6, which 
may imply an increase in the intergalactic medium neutral fraction. From the minimum contribution 
of nebular emission required to fit the observed SEDs, we place an upper limit on the escape fraction 
of ionizing photons to be ~ 0.6 at 2; = 5.7 and ~ 0.9 at z = 6.6. We also compare the stellar 
populations of our LAEs with that of stacked HST/WFC3 z-dropout galaxies. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: evolution — galaxies: 
high-redshift — galaxies: stellar content — 



1. INTRODUCTION 

Lya emitters (LAEs) are a galaxy population com- 
monly seen at high redshift. The high number density 
of LAEs indicates their importance in the evolutionary 
history of galaxies. Since galaxies with strong Lya emis- 
sion can be identified as LAEs even when their contin- 
uum emission is too faint to be detected in broadband 
imaging, LAEs provide an opportunity to probe low-mass 
galaxies with active star-formation. Some of them may 
be building blocks of more evolved galaxies. 

Surveys of LAEs have been made primarily 
with narrow-band inr a ging to isola. t e Lya emission 
(e.g. iHu et al.l [l99l iRhoads et all [2000: live et al.l 
|2006( ). and until now over a thousand LAEs have 
been photometrically s elected a r id /or spectroscop- 
ically identified (e.g., iHu et al.l 120021 : lOuchi et al.l 
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Studying the stellar population of LAEs is essential to 
understand their physical nature and to reveal the rela- 
tionship between LAEs and other h igh-redshift g alaxies 
such as Lyman break galaxies (e.g.. iSteidel et aLl ll996'). 
At z ~ 3 - 5, much progress has been made recently, 
and it has been revealed from large samples with multi- 
band photometry that most LAEs arc small galaxies with 
masses 10® - lO^Mc^ and young ages < 10® yr (e.g. 
Gawiser et all 120071: iNilsson et all 120071: iPirzkal et ahl 
while some ha. v e large stellar masses of ~ 10^° Mr 



,e.g. iLai et"all 120081 : iFinkelstein et al.ll2009l : lOno et all 

mm . 

At higher redshifts, a few studies have reported on stel- 
lar populations of LAEs, but they are all based on a 
very sma l l samp le and consensus has not been reached. 
ILai et all (|2007D have studied stellar populations of three 
bright LAEs at z = 5.7 detected in rest-frame opti- 
cal wavelengths, and found that they have high stellar 
masses of ^ 10^°Mq, as old ages as the Univer se at their 
redshifts, and some amount of dust. lOuchi et"al . (2009^ 
have reported the discovery of a giant LAE at z = 6.595, 
Himiko, with a S pitzer /Infrared Array Camera (IRAC; 
iFazio et al.l 12004 ) counterpart, and estimated its stellar 
mass to be (0.9 - 5.0) xIO^^Mq. It should be noted that 
these four LAEs are all detected in at least one IRAC 
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band, which means that they are not typical LAEs but 
rare, massive-end objects am ong the overa h LAE pop- 
ulation. On the other hand, iPirzkal et alj ()2007l ) have 
studied three faint LAEs at 5.2 < z < 5.8 found by 
Hubble Space Telescope (HST)/ Advanced Camera for 
Survey (ACS) slitless spectroscopy in the Hubble Ul- 
tra Deep Field (HUDF), to show that they are all very 
young (~ several xlO^ yr) with low masses (~ 10^ - 
10^ M(7) ) and smal l dust extinctions {Ay = - 0.1). 
ICharv et al.l ()2005l ) have studied a gravitationally-lensed 
LAE at z = 6.56, HCM6A, and reported that it has a 
stellar mass of ^ 10^ Mq and very young age (~ 5 Myr). 
Although these two studies may be picking up typical 
LAEs at each redshift, it is hard to draw robust conclu- 
sions from such small numbers of o bjects. 

Recentlv. lOuchi et al.l ()2008[ ) and lOuchi et al.l (|2009a[ ) 
have constructed the largest available sample of z = 5.7 
and 6.6 LAEs in an about 1 deg^ of the Subaru/ ZMM- 
Newton Deep Field (SXDF) from deep optical broadband 
and narrowband data. About 65% of the field is also cov- 
ered by deep JHK images taken with the United King- 
dom Infrared Teles cope (UKIRT)/ W ide Field Infrared 
Camera (WFCAM; ICasah et al.l [2001 from UKIRT In- 
frared Dee2_Ski^_Su^e5^KIDS^^ Survey 
(UDS: lLa^ence et al.ll2007HWarren eTaIll2007[ ). and 3.6 
- 8.0/^m images taken with the Spitzer/IRAC from the 
Spitzer legacy survey of the UDS field (SpUDS; Spitzer 
Proposal ID #40021; PI: J. Dunlop). In this paper we 
use these wide-field, multiwaveband survey data to study 
the stellar population of typical LAEs at z — 5.7 and 6.6 
from stacking a large number of faint objects (165 at 
z = 5.7 and 91 at z = 6.6). At each redshift the stacked 
object is detected in several broadbands including the 
3.6/^m band, enabling us to place meaningful constraints 
on stellar population parameters. We also examine the 
stellar population of z-dropout galaxies recently discov- 
ered in the deep HST/ Wide Field Camera 3 (WFC3) 
data of the HUDF, using the stacked spec t ral en ergy dis- 
tribution (SE D) obtained bvlL abbe et al.l (|2010[ ) from 14 
objects of the 
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Oesch et al. (201( 1) sample. 



Schaerer fc de Barroi ([200l . careful SED 



fittings including nebular emission are required to accu- 
rately determine stellar population parameters fo r high- 
z galaxies (see also Zackris son et al. 2008: Raite r et al.l 
[20101: iSchaerer fc de Barrosll20lOir ~ Nebular emission is 
produced in the interstellar medium (ISM) ionized by hot 
stars. Hence, the nebular emission of a galaxy increases 
as the fraction of ionizing photons absorbed by HI gas in 
the ISM increases, or equivalently, as the fraction of ion- 
izing photons escaping from the galaxy, /^°", decreases. 
To see the infiuence of nebular emission on the deter- 
mination of stellar population parameters, we examine 
SED models for two extreme /g°" values: /p™ = 1 where 
the SED is determined solely by stellar emission as most 
previous studies assumed, and /^°" = where the con- 
tribution of nebular emission to the SED is largest. We 
find that the latter models generally give a very good fit 
to the observed SEDs, especially the z = 5.7 one which 
is fit by the latter significantly better. 

In addition, we treat as a free parameter in our 
SED fitting to place rough constraints on /^™, because 
redshifts of z ~ 6 - 7 are close to the end of cosm ic reion- 
ization fe.g.. lFan et al|[200l [Becker et al.|[2007[ ) and the 



of galaxies is a key parameter which determines the 
ionizing photon budget. At redshifts up to z ^ 4, the /, 
of galaxies has been estimated or constrained by a vari- 
ety of methods using, e.g., FUV spect ra and narrow-band 
images of star - formi ng galaxies (e.g., Steidel et al.ll200ll 
iShaplev et all 120061 : llwata et al.ll2009h and the distribu- 
tion of neutral hydrogen column densities i n the after 



glow spectra of long du ration GRBs (e.g., iChen et aD 
120071: iFvnbo et aIll2009D . ahhough the results do not 
agree well with each other. For example, iShaplev et al.l 
(pool) have derived > 0.65 for two z - 3 LBGs with 
detected Lyman continua and J^™ = 0.14 f rom a com- 
posite spectrum of 14 z ~ 3 LBCJ*^. while iChen et al.l 
()2007[) have placed an upper limit of < 0.075 us- 
ing a compiled sample of 28 GRBs (see also lFvnbo et al.l 
[2001 . 

At redshifts as high as z 6 - 7, most constraints 
are based o n the UV luminosit y de nsity of galaxies. 
For ex ample, lOuchi et al.l (|2009bD and iFinkelstein et al.l 
(|2010[) inferred the production rate of ionizing photons 
in galaxies from the observed UV luminosity density of 
z-dropout galaxies and obtained lower limits to to 
keep the intergalactic medium (IGM) ionized. However, 
this method has a number of uncertainties such as the 
UV luminosity functio n of galaxies and the d umpiness 
of the IGM. Recently, iBouwens et al.l (|2010b[ ) proposed 
that the very blue UV color of z-dropout galaxies they 
found may be due to weak nebular emission and hence 
high /o°c , because strong nebular emission makes the UV 
color too red. This is interesting in relating the SED to 
/p™. However, as they already state, is not uniquely 
determined from the UV color but it also depends on the 
age of the stellar population. In this paper we show that 
inclusion of a UV-to-optical color reduces this degener- 
acy and provides constraints on /^°" for our LAEs. 

The outline of this paper is as follows. After de- 
scribing the imaging data used in this study in Section 
2, we produce stacked images in Section 3. Our SED 
fitting method is presented in Section 4. In Section 
5, we present and discuss our SED fitting results and 
constraints on /g°". A summary is given in Section 6. 
Throug hout this paper, we use magnitudes in the AB 
system (jOke fc Gunnlll983D and assume a fiat universe 
with (rj„, f^A, h) = (0.3, 0.7, 0.7). The age of the uni- 
verse is ~ 0.98 Gyr at z = 5.70 and ~ 0.82 Gyr for 
z — 6.56. 

2. DATA 

Deep BVRi' z' images of the SX DF were taken with 
Suprime-Cam (jMivazaki et al.l[2002| ) on the Subaru Tele- 
scope b y the Suhaiu/ X MM- N ew ton Deep Survey pr oject 
(SXDS: iFurusawa et al.l [20081). lOuchi et al.l ([200l and 
lOuchi et al.l ( 2009aD combined this public data set 
with their own imaging data taken with Suprime-Cam 
through two narrowband filters, NB816 (Ac — 8150A, 
FWHM = 120A) and iVS921 (A^ = 9196A, FWHM 
= I32A), and constructed samples of 401 z = 5.7 and 
207 z = 6.6 LAEs over a sky area of ~ 1 deg^ by ap- 
plying the following selection criteria: (i) existence of a 
narrowband excess, (ii) no detection in any bandpasses 

They have assumed the intrinsic UV to Lyman continuum 
flux density ratio to be 3.0, and corrected for IGM absorption. 
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TABLE 1 

Median-Stacked Lyman Alpha Emitters 
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Fig. 1. — Median-stacked images of LAEs at 2 = 5.7 (top) and z = 6.6 (bottom) taken by the Subaru Suprime-Cam {B, V, R, i' , NB, 
z'), the UKIRT WFCAM (J, H, K), and the Spitzcr IRAC (3.6/im, 4.5Atm, 5.8Atm, 8.0/im). The z = 5.7 LAE is detected (brighter than 
the 3(T magnitude) in i' , NB816, z' , and 3.6^m, and nearly detected in J. The z = 6.6 LAE is detected in A'^_B921, z' , and 3.6^m. Each 
panel is 6" X 6" in size, or ~ 35 (33) kpc at z = 5.70 (6.56). 

sizes of the J, H, K images match those of the optical im- 
ages (FWHM K, l."0). The Scr limiting magnitudes over 
a 2"-diameter aperture are calculated to be 24.5, 24.2, 
and 24.4 in the J, H, K bands, respectively. Because the 
zero-point magnitudes for the J, iJ, K images are given 
in the Vega system, we convert them into A B magnitudes 
using the offset values given in Table 7 of iHewett et al.l 
(12001 . 

The SpUDS covers 0.65 deg^ of the overlapping area 
of the SXDS and UDS fields. This 0.65 deg^ area corre- 
sponds to an effective survey volume of 6.0 x 10^ Mpc^ 
for z = 5.7 LAEs and 5.2 x 10^ Mpc^ for z = 6.6 LAEs, 
respectively. All of the SpUDS IRAC images are geomet- 
rically matched to the optical images. We calculate the 
3(7 limiting magnitude over a 3"-diameter aperture to be 
24.8, 24.5, 22.7, and 22.6 in the 3.6, 4.5, 5.8, and 8.0Aim 
IR AC bands, respectively. 

3. STACKING ANALYSIS 

In this paper, we only analyze LAEs in the overlapping 
area of 0.65 deg^ where the Suprime-Cam, WFCAM, and 
IRAC data are all available. We perform IRAC Channel 
1 photometry with a 3"-diameter aperture at the position 
of LAEs in the narrowband images, using the IRAF task 
phot. Among a total of 189 (106) LAEs at z = 5.7 
(6.6), 165 (91) are found to be fainter than the IRAC 
3.6/im-band 3(7 magnitude (i.e., 24.8 mag), which means 
that they are neither rare massive-end objects among 
the overall LAE population, nor confused by neighboring 
objects. We make their median-stacked multi-waveband 
images separately for the two redshifts. Their cutouts 
are shown in Figure [T] 

Among those not used for stacking, some seem to have 
counterparts in the IRAC 3.6/xm image. However, they 
have not been spectroscopically confirmed , except for 
the gia nt LAE Himiko already studied by lOuchi et al.l 
(l2009aD . We will discuss elsewhere stellar populations 
of these bright LAE candidates after we confirm their 
redshifts by spectroscopy. 

3.1. Photometry 

We perform BVRi'z'JHK photometry with a 2"- 
diameter aperture at the position of the LAEs in the 
narrowband images, using the IRAF task phot. We then 
convert them into total magnitudes by subtracting aper- 
ture correction termj^. which are evaluated for bright 

Aperture correction terms [ABmag] are 0.147(-B), 0.115(y), 
0.146(R), 0.202(i'), 0.156(2'). 0.196(J), 0.208(/i'), O.m(K). 



Note. — All magnitudes are total magnitudes. 99.9 mag 
means negative flux densities. Magnitudes in parentheses are 
3cr limiting magnitudes. 

ArB816 for z = 5.7 LAE, NB921 for z = 6.6 LAE. 

Rest-frame UV continuum magnitudes at (observed) effec- 
tive wavelengths of 8303A {z = 5.7) and 9445A (2 = 6.6) after 
correction for Lya emission and the IGM absorption, derived 
from NB8W and i' magnitudes (2 = 5.7) and A'^i3921 and 2' 
magnitudes (2 = 6.6). 

Derived from rricon and J magnitude using eq. |[TJ 

Derived assuming that the redshifted Lyo line is at the cen- 
ter of the NB filter. 

" Corresponding to the central wavelengths of NB816 and 
NB921. 

blueward of the Lyman limit, and (iii) existence of a spec- 
tral break due to IGM absorption. The z = 5.7 LAEs 
have L(Lya) > 3 x lO'^^ g^g g-i and EW(Lya) > 27 A 
(|Ouchi et al.ll2008| ). and the z = 6.6 LAEs have L(Lya) 
> 3 X lO''^ erg s'^ and EW(Lya) > 14 A (|Ouchi et al.l 
[20T0I) . 

About 77% of the SXDF Suprime-Cam field was 
imaged in the J, H, and K bands with the wide- 
field near infrared camera WFCAM on the UKI RT in 
the UKIDSS/UDS project (|Lawrence et aLllMrzl ). The 
UKIDSS/UDS is underway, and we use Data Release 5 
for this study. We align the J, K images with the 
SXDS optical images using common, bright stars, and 
then smooth them with Gaussian filters so that the PSF 
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and isolated point sources in each band. To evaluate the 
aperture correction term for the two narrow bands, we 
measure fluxes for bright point sources in a series of aper- 
tures from 2" up to 6" with an interval of 0.1". Since we 
find that the fluxes level off for > 5" apertures, we define 
the difference in magnitude between the 2" and 5" aper- 
ture magnitudes as the aperture correction term. For the 
Spitzer/IRAC four bands, we measure 3"-diameter aper- 
ture magnitudes for each LAE and converted them to to- 
tal magnitudes by applying the aperture correction given 
by Multiwavelength Survey by Yale-Chile (MUSYC) sur- 
VC5Q. The correction values are 0.52, 0.55, 0.74, and 0.86 
mag for 3.6/im, 4.5^m, 5.8^m, and 8.0/xm, respectively. 
We measure the limiting magnitude for each band by 
making 1000 median-stacked sky noise images, each of 
which is made of (165 for z = 5.7 and 91 for z = 6.6) 
randomly-selected sky noise images. Table [T] summarizes 
their photometry, UV spectral slope /? (sec Section 1X2)) . 
average /(Lya), i(Lya), and EW(Lya). We find that 
the both stacked objects are brighter than the 3a mag- 
nitude in the 3.6/^m band. 

Since our i'-band photometry for the z — 5.7 LAE and 
z'-band photometry for the z = 6.6 LAE are contami- 
nated by Lya emission and IGM absorption, we derive 
the emission-free continuum magnitude mcon at 2± SSOSA 
for z = 5.7 and at ~ 9445A for z = 6.6, respectively. 
We obtain mcon = 27.22 ± 0.14 from NB816- and i'- 
band photometry, and mcon — 27.43 ±0.57 from NB921- 
and z'-band photometry, taking account of the contribu- 
tions of Lya emission and IGM absorpt ion (lMadaulll995D 
to th e photometry in each bandpass ([Shimasaku et ahl 
We use mcon instead of the i'-band (z'-band) 
magnitude for the z = 5.7 (6.6) LAE to derive the UV 
spectral slope [3 in Section 13.21 and constrain the stellar 
population in Section S) 

3.2. UV Spectral Slope 

We estimate the slope of the rest-frame UV continuum, 
/3, from two broad-band magnitudes, mi and m2. 
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mi — m2 
2.51og(Ai/A2) 



- 2. 



(1) 



where Ac and Ac are the central wavelengths of the two 
broadbands. For the z ~ 5.7 LAE, we set (mi,m2) = 
(mcon,"^j) and obtain /3 = —2.9 ± 1.0, where 

^con 

the continuum magnitude at 8303 A (see Table 1) and 
A;! = 12500 A ()Tokunaga et all [200I . These wave- 
lengths correspond to rest-frame 1240 A and 1870 A. 
respectively. Similarly, /3 — —3.0 ± 2.7 is obtained for 
the z = 6.6 LAE from the continuum magnitude at 9445 
A and mj, corresponding to 1250 A and 1650 A, respec- 
tively. Thus both have a very blue slope, although the 
uncertainty is large especially for the z = 6.6 object. 

Figure [3] plots (3 against the rest-frame UV absolute 
magnitude, A/uv> for our objects together with LAEs at 
two lower redshifts (z — 3.1 and 3.7) and LBGs at z ~ 4, 
6, an d 7. For the LA Es we calculate /3 using our own 
data (jOno et al.ll20Tot) : ^ = -2.6 ± 0.2 for z = 3.1 and 
(3 = —2.0 ± 0.2 for z = 3.7 from R and z' magnitudes, 

http: //data. spltzer . caltech.edu/popular/simple/ 
20070601_enhanced/doc/00README_photometry 
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Fig. 2. — UV continuum slope /3 versus rest-frame UV absolute 
magnitude Afyv. The filled circles and open triangles represent 
LAEs and LBGs, respectively, colored according to their redshifts. 
The red and gray circles are our LAEs at ^ = 5.7 and 6.6, respec- 
tively, where the Muv of the 2 = 6.6 LAE has been shifted —0.15 
for clarity. Also shown are the values for stacked LAEs at 2: = 3.1 
(cyan circle) and 2 = 3.7 (blue circle) t Ono et al.ll201lil 'l. 2 ~ 4 
LBGs (bl ue triangles: [B ouwens et al. 201QB), 2 ~ 6 LBGs (red 
triangles: IBouw ens et al.||2009| ), and 2 ~ 7 LBGs (gray triangles: 
IBouwens et al.ll2010bl : ILabbe et al.ll20ia ). 

whose rest-frame wavelengths are 1580 A and 2190 A for 
z = 3.1 and 1390 A and 1930 A for z = 3.7. For the 
LBGs at z '-^ 4 6, and 7, w e take the values given in 
IBouwens et al.l (|2009l . [2010bf ). In addition, we calculate 
/3 = -2.9 ±0.4 for the z ~ 7 stacked LBG fro m its J and 
H magnitude^^ given in ILabbe et all dMl . 

We see a weak correlation in the LBGs of individual 
redshifts that fainter objects have smaller /?, i.e., bluer 
spectra, and a tendency that /? at fixed Muv appears to 
become smaller with redshift. A similar correlation can 
be seen for the LAEs at z = 3.1 and 3.7 if we assume 
that evolution between these redshifts is negligible and 
treat them collectively as z ~ 3.5 objects. In comparison 
with this possible correlation seen in LAEs at z ~ 3.5, 
the LAEs at z = 5.7 and 6.6 are offset toward smaller 
/3, although the significance is at most la levels due to 
the large uncertainties in the /? measurement especially 
at z = 6.6. This offset, if real, might suggest that the 
stellar populations of typical LAEs at z ~ 6 - 7 are 
younger, more metal-poor, and/or with less dust than 
those at z ~ 3.5. 

Moreover, the LAEs at z = 5.7 and 6.6 are slightly 
below the /? - Muv correlation of z ~ 6 and 7 LBGs. 
Although the difference between LAEs and LBGs at 
z ^ 6 — 7 is within la uncertainties, it might im- 
ply that the stellar po pulation is slig h tly di fferent be- 
tween these galaxies. iShaplev et al.l ()2003() have re- 
ported that the UV spectra of z ~ 3 LBGs become 
bluer with increasing Lya equivalent width from large 
negative value s (strong absorptio n) to large positive val- 
ues (see also. iKornei et al.]|2010[ ). Sim ilar trends have 
also b e en found for z ^ 4 LBGs bv iPentericci et al.l 
(l2007l) . iVanzella et all ([2009I) . and IStark et al.l (I2010D . 



The rest-frame wavelengths are 1580 A and 1950 
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This trend might continue to z ~ 6 — 7, since there seems 
to be a hint of systematically bluer continua for LAEs 
than LBGs in Figure [2] 



4. SED FITTING 



(e.g.. 
[2008 



Sin ce LAEs tend to be young star-forming galaxies 



Gawiser et aD 120071 : iPirzkal et aDl2007tlLai et al.l 
Quo et al. f l2010ir it is worthwhile to consider neb- 



ular emission in population synthesis modeling. Here we 
make model SEDs in two extreme cases: = where 
ionizing photons arc totally converted into nebular emis- 
sion, and /^™ — 1 where all ionizing photons escape from 
the galaxy. We call the former the "stellar -|- nebular" 
case, and the latter the "pure stellar" case. We calculate 
nebular spectra (lin es and continua) basica.lly fo llowing 
the manner given in iSchaerer fc de Barrosj (|2009[) . 

We use t he stellar population sy nthesis model of 
GALAXEV (jBruzual fc GharlotJiM)! hereafter BC03) 
to produce st ellar SEP J*^. adopting Salpeter's initial 
mass function (|Salpetedll955[ ) with lower and upper mass 
cutoffs of 0.1 and IOOMq. We assume constant star for- 
mation historjH and consider two stellar metallicities 
Z = 0.2Zq and 0.02^0. 

Nebular emission is calculated under the assumption 
of electron temperature Tc = 10^ K, electron density 
TT-c = 10^ cm""^, and case B recombination. We include 
H recombination lines from Balmcr, Paschen, and Brack- 
ett serieJ*^. W c calculate H/3 line luminosity by (e.g., 
lOsterbrock~Ferland 2006,) 



L(H^) [ergs- 



4.78 X 10" 



-13 , 



1 - f " 



(2) 



where iVLyc is the number of ionizing photons produced 
per second. We do not consider absorption of ioniz- 
ing photons by dust. The luminosities of the other 
H recombination lines are computed from i(H/3) us- 
in g the table of relat i ve int ensities of these lines given 
in iStorev fc Hummed (|1995l ). We also include nebular 
lines from non-hydrogen s using the empirical re lative line 
intens ities compiled by lAnders fc Fritze-v. Alvenslebenl 
(|2003[) , on the assumption that the gaseous metallicity is 
equal to the stellar metallicity. 

Nebular continuu m emission is calculated by (e.g., 
iKrueger et allll995f ) 



(total) 

7^ 

aB 



'1 - f" 



Lyc I 



(3) 



where is the case B recombination coefficient for hy- 
drogen. The continuum emission coefficient --y^'"''''^ ^ con- 
sidering free-free and free-bound emission by H, neutral 

We do not use new population synthesis models which in- 
clude t4iermallypuka giant branch (TP-AGB) stars 



(e.g. . IMarastori'2005': ^ Eruzuai2007l) . because LAEs tend to bo very 
young and the contribution of TP-AGB stars should be negligible 
(e.g.. lUno ct al..,20ia) . 

Since most LAEs are very young, constant star formation 
history (SFH) is a reasonable approximation. Younger ages will 
be obtained if exponentially decaying SFH is assumed, while older 
ages will be obtained for smoothly-risi ng SFH, which h as been 
recently applied by several authors (e.g.. [Stark et aLll2009IV 

^® Wc do not include Lya line, since Lya photons are resonantly 
scattered by neutral hydrogen and its strength is quite uncertain. 
Instead, we derive the Lya-free continuum magnitude rrtcon from 
i'- and A'^B816-band photometry for the z = 5.7 LAE, and z'- 
and A'^i3921-band photometry for the z = 6.6 LAE, as described in 
Section 13.11 



He, and singly ionised He, as well as the two-photon con- 
tinuum of H, is given by 



(total) 



^(HI) , ^(2q) , ^(Hel) »(HC+) (HcII) "(He ) 

^ ^ „(H+) ^ ''^ n(H+) ■ 

(4) 

The emission coefficients 7^,'^ {i = HI, 2q, Hel, Hell) 
for wavelengt hs below an d above l/im are tak en from 
tables 4-9 of lAlled (|1984[ ) and tables I and II of lFerlandl 
(|1980l ). respectiveljO. The abundance ratios are set to 
be rt(He+)/n(H+) = 0.1 and n(Hc++)/n(H+ ) = (e.g., 
iBrown fc MathewsjUlTOl : IKrueger et al][l995l ). 

For the dust extinction of output stellar spec tra, we use 
Calzetti's extinction law (jCalzetti et al.l [20001 ) and vary 
E{B — V)i, as a free parameter over and 1.50 with an 
interval of 0.01. For the dust extinction of nebular emis- 
sion, we assume E{B — V)gas — E{B — F)* as proposed 
bv lErb et al.l (l20oiR . 

Note that we do not consider the effect of dust ex- 
tinction on In other words, we assume that a 
Lyman continuum photon cither ionizes a neutral hy- 
drogen atom or escapes into the intergalactic medium 
through, e.g., holes in the gas. Dust extinction for Ly- 
man continuum emission c an substanti ally reduce the 
hydrogen-ionizing flux (e.g.. Ilnoiiil200ll ). However, the 
exact level of Lyman continuum extinction is difhcult to 
assess even for galaxies in the lo cal Universe, and much 
more so at higher redshifts (e.g.. iZackrisson et al.l[2b08i) . 
although some authors suggest tha t the effect is very 
small for high-redshift galaxies (e.g ., iGnedin et al.ll2008l : 
iRazoumov fc Sommer-Larsenll20Tor ) 

In Section WM we place upper limits to /^°" from the 
minimum luminosity of nebular emission required to re- 
produce the observed SEDs. Those are regarded as con- 
servative upper limits, since a non-zero fraction of Lyman 
continuum photons not converted into nebular emission 
will in practice be absorbed by dust before escaping into 
the IGM. 

We perform the sta ndard SEP fittin g method (for de- 
tails, see Section 3 of lOno et aT]l2010( ). We make a large 
set of stellar-mass-normalized model SEPs, varying age 
and dust extinction. We then redshift them to z = 5.70 
and 6.56 and convolve them with bandpasses to calculate 
flux densities. For each object, wc search for the best- 
fltting SED that minimizes separately for /^°" = 
and 1 and separately for Z = 0.2Zq and O.O2Z0. Since 
stellar mass Mgtar is the amplitude of a model SEP, we 
obtain the best-fltting Mgtar by solving dx^ /dMgtar = 0. 
The errors in the best-fit SEP parameters correspond to 
the la confidence interval (Ax^ < 1) for each parameter. 

Since our i'-band flux density for the z = 5.7 LAE and 
z'-band flux density for the z — 6.6 LAE are contami- 
nated from Lya emission and IGM absorption, we use 
'Ticon calculated in Section 13.11 instead of these flux den- 
sities. We do not use short wavebands (i.e., BVR for 
z = 5.7 and BVRi' for z ~ 6.6), since they suffer from 
the IGM absorption shortward of the Lya wavelength. 



1''' We assume 7^'^'^^ = 
ISchaerer fc Vacca|[l99g ). 



and 7, 



(Hoi) 



(HI) 



at A > l/im (e.g 
E(B-V)^,/0A4:, which 



We find t hat adopting E{B-V )sb. 
is proposed by ICalzetti et al.l I I2000I) , does not significantly change 
our results. 
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the amount of which considerably differs among the lines 
of sight. Thus, model SEDs arc fitted to the observed flux 
densities in mcon, -z', J , H, K, 3.6/im, A.Bfim, 5.8/im, and 
8.0/im-band for z = 5.7, and to the same bands except z' 
for z = 6.6. The free parameters in the fitting arc stellar 
mass, age, and dust extinction. The degrees of freedom 
are six for z = 5.7 and five for z = 6.6. 

5. RESULTS AND DISCUSSION 

5.1. Stellar Populations of Lyman Alpha Emitters 
at z ^ 6 -7 

Table [2] summarizes the results of the SED fitting for 
our LAEs. First, wc find that in both f^f^ = and = 
1 the best-fit models for Z ~ 0.2Zq and 0.02^0 are 
very similar to each other. This implies that the changes 
in the stellar and nebular emission spectra over 0.02 < 
Z/Zq < 0.2 are not large enough to significantly alter 
the best-fit parameters. In what follows we concentrate 
on the results for Z = 0.2Zq for simplicity. 

Next, we find that the best-fit models are extremely 
different between = and 1. For = the best- 
fit models have relatively small stellar masses, 3 x IO^Mq 
(z = 5.7) and 1 x 10® A/0 {z ~ 6.6), and young ages, 3 
Myr (z = 5.7) and 1 Myr {z = 6.6). On the other hand, 
for = 1 the stellar masses are more than one order 
of magnitude higher, 5 x 10^ Mq and 3 x 10^^/0, and 
the ages arc more than two orders of magnitude older, 
300 Myr and 800 Myr. As we see below, these great 
differences are related to how to account for the observed 
bright IRAC magnitudes, or equivalently the red UV-to- 
optical colors. In contrast to the stellar mass and age, 
the dust extinction is consistent with E{B — F)* = in 
both cases for both objects, suggesting that typical LAEs 
at z 6 - 7 are nearly free from dust extinction. 

Figure |3] compares the best- fit model spectra with the 
observed flux densities. The bottom panels show the 
results for /g™ = 1, the "pure stellar" case. For z = 5.7, 
the best-fit model matches the observation at Aobs ^ 2fim 
but undershoots the data point at 3.6/im. This offset 
will be reduced if older ages or larger E{B — T^)^ values 
are adopted, but such models will then not fit the blue 
UV spectra. The observed SED of the z = 6.6 LAE 
resembles that of the z = 5.7 LAE, with blue UV colors 
and detection at 3.6//m. In addition, the z = 6.6 LAE has 
a relatively bright flux density at 4.5/im as well although 
less than 3cr detection. The overall shape of the z = 6.6 
SED is reproduced well by a pure stellar model. 

The top panels of Figure [3] show the results for /^™ — 
0, i.e., the "stellar -|- nebular" case. We find that for 
z — 5.7 the discrepancy at 3.6/im seen in the bottom 
panel almost disappears thanks to strong nebular emis- 
sion lines such as [Oiii] and II/3 contributing to this 
bandpass; the stellar emission has only a minor contri- 
bution. Thus, adopting a very young stellar population 
with strong nebular emission can simultaneously fit the 
observed blue UV color and the red UV-to-optical color. 
For z = 6.6 as well, the best-fit model with /^™ — 
reproduces the observed SED well in a similar manner. 

Then, which of — 1 and = is more favored? 
For z = 5.7 it is easy to answer this question. As shown 
in Table [H the best-fit for f°sc = 0, 7.75, is signifi- 
cantly smaller than that for /g™ = 1, 10.3. Indeed, the 
/osc ~ model fits the 3.6/tm data point and blue 3.6/im 



- 4.5/im color far better than the /g™ = 1 model. In ad- 
dition, the small photometric errors in the mcon, z' , and 
J magnitudes do not permit old-age, high-mass models 
which fit the 3.6/im data but instead give red UV colors. 

For z = 6.6, on the other hand, the best-fit /g™ = 1 
model gives almost the same (2-9) as the /g°" = 
model (2.5), and we cannot conclude which is favored 
solely from the values. This is partly because the 
larger photometric errors in the J and shorter band- 
passes, and flatter 3.6/im - 4.5/im color permit old, mas- 
sive models with red UV colors. However, typical LAEs 
at z 3 - 5 have been consistently found to be young and 
low- mass objects (e.K. , | Gawiscr et al. 2007t iPTrzkal et abl 
[2OO7I: iLai et al.l [2001 fOno et all I2O1O0 . and our study 
shows that this trend continues at z = 5.7. Combining 
this fact with an argument that it seems unnatural that 
LAEs at earlier epochs are older and more massive, we 
take the model with = for the z = 6.6 LAEQ- 

We conclude that typical LAEs at z ~ 6 - 7 have 
low stellar masses of (3 — 10) x lO^Af0, very young 
ages of 1 — 3 Myr, and negligible dust extinction. 
Wc thus propose that they are candidates of galaxy 
building blocks at an early stage of galaxy formation. 
This proposal can be regarded as an extension toward 
higher redshift of a similar idea which has been pro- 
posed for low- red shift LAEs based on apparent mag- 
nitu des and sizes (iPascarelle et al.l 119961 for z ~ 2; see 
also lOuchi et al.l [200f who found that z ~ 5 LAEs are 
UV f aint) and stellar population analysis (jGawiser et al.l 
120071 for z ^ 3). This proposal is also consistent with 
low dark-halo masses estimated for LAEs from clustcr- 
ing a nalysis (|Gaw iser et al. 2007 for z ^ 3. iGua ita et al] 
[2OIOI for z ^ 27 1buchi et7alll2010l for z - 3 - 7). 

5.2. Constraints on Lya Escape Fraction 
We estimate for our LAEs the Lya escape fraction, 

./esc ' 

rLya ^ -^obs(LyQ) 

/esc T [T \ ? \^) 

Lint (Lya) 

where Lobs (Lya) is the observed Lya luminosity and 
Lint (Lya) is the intrinsic Lya luminosity computed from 
the SFR on the assumption of case B using Lint(Lya) 
[ergs-^] = 1.1x 10^^ SFR [M0 yr"!] (|Brocklehurstill97ll: 
IKennicut3ll998D . As is easily noticed, /^^" corresponds 
to the fraction of Lya photons produced in the galaxy 
which escape from absorption by the galaxy's ISM and 
absorption by the IGM at the galaxy's redshift. Because 
the dust extinction of our LAEs is negligibly small (Sec- 
tion l5.ip , Lya photons propagating in the ISM are mostly 
just scattered by HI gas without absorption. In this case, 
/^^" depends mostly on the strength of the IGM absorp- 
tion. 

Substituting the observed Lya luminosity and SFR de- 
rived from our SED fitting, we obtain /^y" ~ 0.36+0:35 
for z = 5.7 LAEs, and /^y" ~ 0.040to^38 for z = 6.6 

We cannot completely rule out the possibility that due to some 
selection effect unique to the z = 6.6 sample, our z = 6.6 LAEs 
are biased toward very massive objects. However, such selection 
effects are very unlikely, since the z = 6.6 sample has been selected 
in a similar manner to the z = 5.7 sample from the same data set, 
and the stacked objects at these two redshifts have similar L(Lya) 
and Afuv, as found in Table [T] 
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Fig. 3. — Best-fit SEDs (curves) and observed magnitudes (filled squares: used for SED fitting, open squares: not used) for the z = 5.7 
LAE (left panels) and the z = 6.6 LAE (right panels). The top panels arc for 'stellar + nebular' (/^°" = 0) models. The red solid curves 
indicate the best-fit SEDs that are the sum of a stellar SED (red dotted curve) and a nebular SED (red dashed curve) . The crosses indicate 
synthesized flux densities in individual bandpasses. The bottom panels are for 'pure stellar' (/c™ = 1) models. 



LAEs, as plotted in Figure [5] We fi nd a tentat i ve de - 
crease in f^l" from z = 5.7 to 6.6. lOno et all ([20lol) 
ha ve found that LAE s at z - 3 - 4 have f^/f ~ 0.1 - 
1 . iHaves et al.l (j20100 have obtained the median value of 
flZ" for z = 2.2 LAEs to be higher than 0.32. These re- 
sults might suggest that the Lya escape fraction of LAEs 
is nearly constant up to 5.7, and then decreases toward 
z = 6.6. This decrease, if real, could be due to an in- 
crease with redshift in the hydrogen neutral fraction of 
the IGM, thereby Lya photons are scattered more fre- 
quently. Indeed, a possible increase in the neutral frac- 
tion has been proposed from an observed decline in the 
nu mber density of bright L AEs from z = 5.7 to 2 = 6.6 
bv lKashikawa et al.l ()2006[ ). 

Further discussion on the neutral fraction is, however, 
difficult given the quality of the fl'/^" measurements. 
Note also that there may be some other mechanisms 



which change fc/c"- For example, unlike our assump- 
tion, Lya emissivity at a given SFR might decrease from 
z = 5.7 to 6.6. Although Figure [S] shows potential use- 
fulness of fc/c") as a means to evaluate the IGM neutral 
fraction near the reionization epoch, much deeper multi- 
band photometry and more detailed SED modeling will 
be needed to obtain a reliable constraint. 

5.3. Comparison with LBGs at z 7 

We fit model SEDs to a z ~ 7 z-dropout composed of 
14 objects recently discovered by the IIST/WFC3 survey 
(jOesch et al.l 120101: iLabbe et al.il2010D . to compare with 
our LAEs. They have been selected from an extremely 
deep (5(7 « 28.6 — 28.7 over 0."25-radius apertures in 
Y105, Ji25; and -ffieo) area of 4.7 arcmin'^ using zgso ~ 
Y105 and Y105 — J125 colors and have an expected redshift 
distribution z ~ 6.5 — 7.5. We consider the observed 
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TABLE 2 

SED Fitting Results for the LAEs at z ■ 



5.7 AND 6.6 



model 


Z 

\Zq\ 


log Mstar 


E(B - y)* 
[mag] 


log{Age) 


log(SFR) 


log(SSFR) 




2 = 5.7 LAE 


stellar + nebular 
pure stellar 

stellar + nebular 
pure stellar 


0.2 
0.2 
0.02 
0.02 


7 4Q+0-20 
7 49+0-37 


00+" "^ 
00+°"^ 
00+° °5 
00+° '"' 


6.50t?ig 
8.46tO:|5 
6.62t?-52 

a +0.25 


qq+l *« 
u.yy_o.46 

o.sit":^^ 

o.87tJ;i 
0.26tO;g5 


-6.50ll;« 

-»-3»_0.48 

o fi-1+1-13 
-8.bi_0 24 


7.75 
10.3 
9.00 
10.4 



z = 6.6 LAE 



stellar + nebular 
pure stellar 

stellar + nebular 
pure stellar 



0.2 
0.2 
0.02 
0.02 



7 05+1-66 
'■^^-0.30 
q 43+0.28 
^■^•^-0.33 

9.57«:i 



'-'■11-0.11 ''■y'^-0.85 

n nq+0-19 s 01+0 00 
0.14+0-?? 8.91+0-«« 



r+2.96 



-0.11 



-0.85 



2.00 
0.63 
1.69 
0.77 



+1.21 
1.66 

+0.71 
0.24 
+ 1.63 
1.27 
+0.74 
0.28 



-5.95 
-8.81 



0.85 
2.86 
+2.33 



i.38_2 42 

sn+2-26 

'•*^-0.00 



2.53 
2.86 



0.00 
+ 1.28 2.46 

3.05 



TABLE 3 

SED Fitting Results for the 2 ^ 



7 2-DROPOUT GALAXY 



model 


Z 

[Ze] 


log Matar 

[Mq] 


EiB - V)* 
[mag] 


log{Age) 

[yr] 


log(SFR) 

[Mo yr-i] 


log(SSFR) 

[yr-^] 




z ~ 7 2-dropout 


stellar -|- nebular 
pure stellar 


0.2 
0.2 


7 79+0.10 
' ■ '^-0.18 
og"l~0.22 

«-*''^-0.30 

8.98«:- 


^■'J^-0.04 

'J-'^^-o.oo 
^■'J'-'-o.oo 


6.16«;^^ 
8.76t°i« 


1 56+^-^'^ 
l-^°-0.68 
O1+0.05 

'^■■'i-o.oi 

r, 09 + 0.05 

"■■'^-0.00 


-6 16+^°'' 

c-r, + 0.67 

o-a^_0.24 


7.35 

6.90 
6.26 



Note. — For the 2 ~ 7 2-dropout, shown are two equally well fitted /^°" 
young model and a very old model. See text for details. 



models, a very 



flux densities in the -Zssoi ^1 



105; 



K , 3.6/im , 



4.5/im band taken from Table 1 of lLabbI' eFall (l2010h . 
assuming a redshi ft of z = 6.88, which is derived by 
iLabbe et al.l ()2010D . This stacked object is similarly faint 
in the rest UV continuum (J125 = 29.6) to our LAEs, 
and thus suitable for compariso n. We do not use tw o z- 
dropout candidates reported by ICapak et al.l (120091) an d 
11 ^ssoLP-dropouts discovered bv lGonzalez et al.l (j2010D . 
because they are too bright (the former have J 23 and 
the latter Jnow ~ 26 — 27.5). In our SED fitting, we 
rule out the Z775-band data although ILabbe et al.l (120101 ) 
used it for their SED fitting, since it suffers from the 
IGM absorption shortward of Lya wavelength, and the 
amount of the absorption differs with the line of sight. 
The free parameters are thus stellar mass, age, and dust 
extinction, and the degrees of freedom arc four. 

This object has also been SED-fitted by 
ISchaerer fc de BarrosI (|2010t) . using stellar popula- 
tion synthesis models with nebular emission. The major 
difference between their and o ur modelings are the 
star f ormation history assumed. ISchaerer fc de Barroj 
(|201Cl[ ) have assumed an exponentially declining star 
formation rate with its e-folding time as a free param- 
eter, while we assume a constant star formation rate. 
For a fair comparison with our results on LAEs, we 
reanalyze the z-dropout SED with our models. The 
best-fit parameters are given in Tabled 

The bottom panel of Figure [S] shows the result for 
/_;™ = 1. We find that the observed SED is well ex- 
plained by a low-mass (log (Mstar[-^o]) = 8.98tQ;27) ^'^'^ 



moderately aged (log (Age[yr]) = 8.76_g3g) model with 
little dust extinction {E{B - V)^ = 0.00+g : g^g|. This 
result is consistent with that of ILabbe et al.l (j2010t) who 
examined /^°" = 1 models alone. 

The top panel of Figure [H] shows the result for /^°" = 0. 
Wc obtain almost the same parameters as the = 1 
model: log (Af^tarlMg]) = 8.83t°:22, log(Age[yr]) = 
8.6lt°j|3, and E{B ~ V)^ = 0.00l°°J. However, we 
find that the data are also well reproduced by a low- 
mass (log(M,tar[M0]) = 7.72t°;;[), extremely young 



(log (Age[yr]) 



6 16+0-30 

1+0.02 



), and almost extinction free 



{E{B - V)^ = G.04lo-u^) model, which are also shown 
in Table [2l Although these two models are extremely 
different, they have almost the same values, and we 
cannot conclude which is more favored. In summary, the 
z-dropout galaxy is either a very young, low-mass ob- 
ject, or a very old, massive object. To be interesting, 
just like in the case of the z = 6.6 LAEs, models with 
intermediate ages and masses are not favored. 

If we take account of the difference in the assumed 
star formatio n history, our results a. r e broa dly consistent 
with those of ISchaerer fc de BarrosI ()2010D . As shown in 
Figure C.2 of their paper, they have derived two equally- 
fitted solutions: very young (^ 4 Myr) and old (^ 700 

^0 The age suggests that the formation redshift of this object is 

> 9.5. 

The age suggests that the formation redshift of this object is 

> 8.5. 
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Fig. 4. — of the best-fit model as a function of age (top: 
2 = 5.7 LAE, bottom: z = 6.6 LAE). The red curves are for 
= and the black curves for = 1. 
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Fig. 5. — Lya escape fraction as a function of redshift. The cir- 
cles (squares) correspond to z = 5.7 (6.6). The filled and open sym- 
bols are for = and 1, respectively, where the open symbols 
have been shifted by —0.1 along the x-axis for clarity. Although 
/esc = 1 is physically unrealistic because Lya photons cannot be 
produced, wc plot the results to show how f^sc°' varies with 
in our method. 



Fig. 6. — Same as FigureO but for the 2 ~ 7 2-dropout galaxy 
IILabbe et al.l20lOl) . In the top panel, two equally well fitted /^°" = 
models are plotted: the 'young' model in red and the 'old' model 
in blue. 
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Fig. 8. — Top panel: age versus stellar mass Mstar- The circles, 
squares, and triangles are for the z = 5. 7 LAE, z = 6.6 LAE, and 
2 ~ 7 dropout galaxy l|Labbe et al.l20Toh , respectively. The results 



for r° 



■ are colored in red, while those for 



1 are shown 



in open symbols. For the z ~ 7 dropout galaxy, we also plot the old 



■I es 



model in blue. Bottom panel: specific star formation 



rate (SSFR = SFR/Mstar) versus stellar mass Mstar- The symbols 
are the same as in the top panel. 

Myr). Note that their values do not match with ours. 
This would be because they have parameterized the star- 
formation e-folding time and metallicity, which we fix, as 
well, and they have considered -8435, Veoej and 1775 data, 
which we do not consider. 

The top panel of Figure [5] plots the best-fit age against 
the best-fit Afgtar for our LAEs and the z ~ 7 z-dropout. 
The best-fit solutions for our z = 5.7 and 6.6 LAEs (red 
filled circle and square) suggest that they are low-mass 
and very young star-forming galaxies. As for the z ^ 7 
z-dropout, both of the two extreme solutions (red and 
blue filled triangles) are plotted. If the younger solution 
is true, the z-dropout is as old and massive as our LAEs. 
On the other hand, if the older solution is true, the z- 
dropout is much older and more massive than our LAEs. 

The bottom panel of Figure [5] shows the specific SFR 
(SSFR) as a function of Mgtar- Our z = 5.7 and 6.6 
LAEs have low stellar masses and very high SSFRs. On 
the other hand, the z-dropout has either almost the same 
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Fig. 9. — X versus the Lyman-continuum escape fraction 
for the 2 = 5.7 LAE (top panel) and the 2 = 6.6 LAE (bottom 
panel). The dashed lines correspond to Xmin ~^ ^ ■ 
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Mstar and SSFR as the LAEs. or far more massive and 
lower SSFR than the LAEs. 

5.4. Constraints on Ly Continuum Escape Fraction 

In the previous sections we consider two extreme val- 
ues for the Ly continuum escape fraction: /g™ = 1 (pure 
steUar) and /g°" = (stellar -|- nebular). Indeed, our 
data are not deep enough to place as a strong constraint 
on as on the other parameters. However, since 
is a very important quantity which controls cosmic rcion- 
ization, in this subsection we perform SED fitting with 
/^°" as an additional free parameter, to try to obtain 
rough constraints on /^™. We vary over and 1 
with an interval of 0.1. In the fitting of the z = 6.6 LAE, 
we search for the best-fit age in the range of < 20 Myr 
following the argument in Section l5.ll that the z = 6.6 
LAE is likely to be very young. 

Figure [9] shows as a function of /p™ for our LAEs. 
For z = 5.7, is nearly constant up to /^°" = 0.5 and 
then starts to increase. Beyond /^°" ^ 0.6 it exceeds 
Xmin + 1' the 1(7 confidence level. We can thus place 
an upper limit of ~ 0.6. For z = 6.6, changes 
little up to as high as = 0.8, and exceeds Xmin + 1 
at around 0.9. Thus is only loosely constrained as 
/esc ^ We also do a similar analysis to the z-dropout 
but obtain no meaningful constraint, since x^ does not 
change larger than unity over the whole /^°" range. This 
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TABLE 4 

SED Fitting Results for the LAEs at 2 = 5.7 and 6.( 
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Fig. 10. — Left panel: comparison of the observed 2' — [3.6] versus z' — J of the z = 5.7 LAE (shaded region: la range; hatched region: 
2(T range) with model tracks from ~ 10® yr to 9 X 10* yr with Z = O.2Z0 and a constant star formation rate, for /g™ = (solid curve 
with red circles corresponding to four ages: 10®, lO'^, 10*^, 9 X 10* yr from left to right), 0.5 (dotted curve with magenta squares), and 1 
(dashed curve with blue triangles). The open symbols represent ~ 10® yr models with = 0.1 — 0.4 (red open ci rcles) and f^g° = 0. 6 
- 0.9 (magenta open squares) from top to bottom. The arrow shows the effects of dust extinction on model colors llCalzetti et al]|2000l '). 

[3.6] versus J125 — Hiao of the z ~ 7 z-dropout IjLabbe et al.ll2010l ) (shaded region: la 



Right panel: comparison of the observed J125 
range; hatched region: 2a range) with model tracks from 



10 yr to 7 X 10 yr with Z = 0.2Zq and a constant star formation rate, for 



J es 



(solid curve with red circles corresponding to four ages: ~ 10®, lO'', 10®, 7 X 10* yr from left to right), 0.5 (dotted curve with 



mage nta squares), and 1 (dashed curve with blue triangles). The arrow shows the effects of dust extinction on model colors l lCalzetti et al.l 
[2000f) . 



is probably because the z-dropout is fit well by an old-age 
model as well (see Section [53]) in which nebular emission 
is so weak that is insensitive to /g™. If we limit the 
age of the z-dropout to < 20 IVIyr, as in the case of the 
z = 6.6 LAE, then we obtain an upper limit of /^°" ^ 0.2. 

Although our LAE data cannot strongly constrain /g™, 
it would be worth presenting the best-fit results for 
= 0.2, at which for the z = 5.7 LAE reaches 
its minimum. As found in Table |4l the best-fit parame- 
ters for = 0.2 overlap well with those for = 
(Table ^ within the errors, indicating that our conclu- 
sions on the stellar populations obtained in Section 5.1 
are relatively robust against the uncertainty in 

The reason for /^°" having an upper limit is that with 
too large /^°", young stellar populations favored by the 



observed blue UV continua cannot produce nebular emis- 
sion strong enough to account for the observed red UV- 
to-optical color. In this sense, the UV-to-optical color is 
critical to constrain /g™. We explain this situation us- 
ing Figure 1101 In this figure, three model tracks for the 
z = 5.7 LAE over ages of 10^ to 9 x 10® yr are plotted 
according to three /^°" values, 0,0.5, 1, in the z' — [3.6] 
versus z' — J plane. We take z' — J as a representative 
of the rest-frame UV color and z' — [3.6] as a rest-frame 
UV-to-optical color bracketing 4000 A. The shaded and 
hatched regions are, respectively, the observationally per- 
mitted 1(T and 2(7 ranges. The z' — J color at a fixed age 
becomes bluer with but the change is modest. In 
contrast, z' — [3.6] increases with /^°" very sensitively for 
young ages. For example, for an age of 10^ yr, which is 



12 



Ono ct al. 



close to the best-fit age of tlie z = 5.7 LAE, z' ~ [3.6] 
becomes rapidly redder with (magenta squares) and 
the model goes out of the hatched region at /^°" ~ 0.3. 
Although the goodness of models should be measured us- 
ing the all magnitude data, this figure demonstrates the 
importance of a UV-to-opti c al color in constraining 

Recently, iBouwens et al.l ()2010b( ) proposed that the 
very blue UV color of z-dropout galaxies they found may 
be due to weak nebular emission and hence high 
because strong nebular emission makes the UV color too 
red. However, as they already state, /g™ is not uniquely 
determined from the UV color but it also depends on 
the age of the stellar population. Inclusion of a UV-to- 
optical color greatly reduces this age-/^°" degeneracy, as 
illustrated in Figure [TOl Basically, with very accurate 
measurements of a UV color and a UV-to-optical color, 
we can obtain a stringent constraint on (for a fixed 
metallicity and IMF). For instance, if /3 of the z-dropout 
was found to be —3 with a great accuracy, low /^°" mod- 
els such as ~ would be ruled out. On the other 
hand, if /? and J125 — [3.6] of the z-dropout were found 
to be —2.8 and respectively, with great accuracies, low 
/g™ models would be favored. 

Several studies have recently constrained lower limits 
of the escape fraction based on the UV lurn inosity density 
of 2;-dropout galaxies: lOuchi et al.l ()2009b[) have obtained 
/e°c ^ 0.2, and iFinkelstein eTall (pOlOl ) have obtained 
/osc ^ The upper limits obtained above complement 
these lower limits, thus narrowing the range permitted. 

Several theoretical studies have argued that could 
vary depending on their host dark halo m asses. How- 
ever, t here is no consensus on this tendency; iWise &: CenI 
(|2009| ) have report ed that /j™ decreases as the ha l o mass 
decreases, while iRazoumov fc Sommer-LarsenI (l201Clfl 
have shown an opposite tendencv (see also. lYaiima et al.l 
I2OIOD . This controversy may be resolved if our method 
is applied to objects in a wide mass range. 

5.5. Contribution of LAEs to the Stellar Mass Density 
and the Cosmic Star Formation Rate Density 

We estimate the contribution from LAEs to the cosmic 
stellar mass density by dividing the total stellar mass of 
LAEs by the comoving volume searched by each narrow 
band. The total stellar mass at each redshift is defined 
as the stellar mass of the stacked LAEs multiplied by 
their number. Assuming /^°" ~ 0, we obtain ~ 8.5 x 10^ 
[Mq Mpc-3] for z = 5.7 LAEs and ~ 1.6 x lO'^ [M© 
Mpc~'^] for z = 6.6 LAEs. These values should be taken 
as lower limits, since we have excluded objects detected 
and/or significantly confused by neighboring objects in 
t he 3.6/xm irnage. 

lEvles et "all (|2007[ ) have summed up the stellar masses 
of the «-drop galaxies they detected, to obtain a lower 
limit to the stellar ma ss density of 2: ^ 6 galaxies of 
- 3 X 10^ [Afo Mpc-3]. iStark et al.l ([2001 have derived 
the mass function of i'-dropout galaxies and obtained a 
stellar mass density of ~ 4.9 x 10^ [Mq Mpc^'^] by in- 
tegrating the mass function brightward of Misoo = —20. 
The stellar mass density of our z = 5.7 LAEs is only 
^ 0.2 — 0.3% of these values. 

For z ~ 7 galaxies, iLabbe et al.l ()2010f ) have obtained 
a lower limit of ~ 3.7 x 10^ [Mq Mpc~^] by multi- 
plying the UV luminosity density of z-dropout galaxies 



given in IBouwens et al.l (|2010aD by an estimated mass- 
to-luminosity ratio. Again, the stellar mass density of 
our z — 6.6 LAEs is as low as ^ 0.4% of this valu e. 

Using SED fitting similar to ours, lOno et al.l ()2010l ) 
have obtained 1.4 x lO'^ [Mq Mpc'^] for z = 3.1 LAEs, 
and 5.2 x 10^ [Mq Mpc^^] for z = 3.7 LAEs. These 
values are about 10—60 times larger than those of z = 5.7 
and 6.6 LAEs, although they assumed /^°" = 1. 

Similarly, we estimate the contribution from LAEs to 
the cosmic star formation rate density (SFRD), where 
the total star formation rate is computed as the star for- 
mation rate of the stacked objects multiplied by their 
number. We obtain ~ 2.7 x 10"^ [Mq yr~i Mpc~^] for 
z = 5.7 LAEs and ~ 1.8 x 10"^ [Mq yr"! Mpc'^] for 
z = 6.6 LAEs. These values are also lower limits to the 
t otal stellar mass density of LAEs. 

IBouwens et al.l (|2007D have estimated the SFRD of 
z ^ 6 galaxies to be ~ 7.2 x 10"^ [Mq yr"^ Mpc"^] 
by integrating the UV luminosit y functi on of i-dropout 
galaxies. Similarly, lOuchi et al.l ()2009b[ ) have obtained 
^ 7.5 X 10""^ [Mq yr~^ Mpc~^] for z ^ 7 galaxies from 
the UV luminosity functi on of z-dropout galaxies (see 
also. lBouwens et alll2010al) . These two values are compa- 
rable to the lower limits obtained for our LAEs at z = 5.7 
and 6.6, suggesting that LAEs are m ajor sources of the 
cosmic star formation at z ^ 6 - 7. lOno et al.l (f2010h 
have obtained ~ 8.2 x lO'^ [Mq yr'^ Mpe'^] for z = 3.1 
LAEs, and ~ 1.3 x 10"^ [Mq yr^^ Mpc^^j foj. ^ ^ 3 7 
LAEs. These values are about 0.5 — 50 times larger than 
those of z = 5.7 and 6.6 LAEs. 

6. CONCLUSIONS 

In this paper, we investigated the stellar populations of 
LAEs at z = 5.7 and 6.6 found in 0.65 deg^ of the SXDF, 
based on deep rest-frame UV-to-optical photometry ob- 
tained from the three surveys: SXDS, UKIDSS/UDS, 
and SpUDS. We made composite images from 165 and 91 
LAEs at z = 5.7 and 6.6, respectively, which are fainter 
than the 3a magnitude in the IRAC 3.6/im band, and 
derived typical SEDs of z ~ 6 — 7 LAEs for the first 
time. We found that their UV continua are as blue as 
those of dropout galaxies at similar redshifts, with UV 
spectral slopes /3 ^ —3, albeit with large photometric un- 
certainties. Fitting stellar population synthesis models 
with and without nebular emission, which is parameter- 
ized by to the multiband data of the stacked objects 
at z = 5.7 and 6.6, we derived their stellar masses, ages, 
and dust extinction. 

Our main results are as follows: 

(i) We find that the stacked LAEs at both redshifts 
are fitted well by = models. The best- 

fit /osc ~ models have low stellar masses of 
(3 — 10) X 10 ''Mq, very young ages of ~ 1 — 3 
Myr, and negligible dust extinction. In these mod- 
els, young stellar populations reproduce the ob- 
served blue UV continua, and strong nebular emis- 
sion redward of 4000 A makes the UV-to-optical 
color as red as observed. 

While we find that the z = 6.6 LAE is also fitted 
similarly well by an old, massive model without 
nebular emission, we do not take this models as 
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the best-fit model, since typical LAEs up to z ~ 6, 
including our z = 5.7 LAEs. have been consistently 
found to be very young and low-mass galaxies. We 
propose that typical z 6 — 7 LAEs arc candi- 
dates of galaxy building blocks at the early stage 
of galaxy formation. 

(ii) We estimate the Lya escape fraction to be f^J^" ~ 
0.36 for the z = 5.7 LAE and f^Z" ~ 0.04 for the 
z = 6.6 LAE with large errors. This decrease from 
z = 5.7 to 6.6, if real, might be due to an increase 
in the neutral fraction of the IGM. 

(iii) We also apply SED fitting to a stacked object from 
z ~ 7 z-dropout galaxy candidates found by a re- 
cent HST/WFC3 deep survey, and find that it is 
either a young, low-mass galaxy with strong nebu- 
lar emission similar to the z = 5.7 and 6.6 LAEs, 



or a very old and massive galaxy with little nebular 
emission. 

(iv) From the constraints on nebular emission models, 
we estimate the upper limit of the Ly-continuum 
escape fraction to be /^°" ^ 0.6 for the z = 5.7 
LAE and ^ 0.9 for the z = 6.6 LAE. We apply this 
technique to the z-dropout galaxies, but obtain no 
meaningful constraints on 
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